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The decomposition of matrices associated to two-qubit and three-qubit orthogonal gates is 
studied, and based on the decomposition the synthesis of these gates is investigated. The optimal 
synthesis of general two-qubit orthogonal gate is obtained. For two-qubit unimodular orthogonal 
gate, it requires at most 2 CNOT gates and 6 one-qubit R y gates. For the general three-qubit 
unimodular orthogonal gate, it can be synthesized by 16 CNOT gates and 36 one-qubit R y and R z 
gates in the worst case. 
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I. INTRODUCTION 



In quantum computing, the algorithms are commonly described by the quantum circuit model |l|. The building 
blocks of quantum circuits are quantum gates, i.e., unitary transformations acting on a set of qubits. In 1995, Barenco 
et al showed that any qubit quantum circuit can be decomposed into a sequence of one-qubit gates and CNOT gates 
Q. The process of constructing quantum circuits by these elementary gates is called synthesis by some authors. The 
complexity of quantum circuit can be measured in terms of the number of CNOT and one-qubit elementary gates 
required. Achieving gate arrays of less complexity is crucial not only because it reduces resource, but it also reduces 
& , errors. 

Decomposition of matrix plays very important role to synthesize and optimize quantum gates. Based on Cartan 
decomposition [3,|j|, the synthesis, optimization and "small circuit" structure of two-qubit gate are well solved [S-Q- 
To implement the general two-qubit gate, it requires at most 3 CNOT gates and 15 elementary one-qubit gates from 
the family {R y ,R z } 0,||. 

Unfortunately, the aforementioned optimal synthesis of most general two-qubit quantum gates have not yet led to 
similarly tight results for three-qubit gates. Based on one of Cartan decompositions for multi-qubit system, Khancja- 
Glaser decomposition (KGD) J4J, Vatan and Williams get the result that a general three-qubit quantum gate can be 
L~ synthesized using at most 40 CNOT gates and 98 one-qubit R y and R z gates (lOJ . Using the modified KGD, the 

results have been improved in [1J|, that is it requires at most 26 CNOT gates and 73 one-qubit R y and R z gates. 
Now the best known result is based on quantum Shannon decomposition (QSD) 12] proposed by Shende, Bullock and 
Markov, it requires at most 20 CNOT gates. According to the result of multi-qubit case, the best known theoretical 
lower bound on CNOT gate cost for general three-qubit gates is 14 [8j. However, no circuit construction yielding 
these numbers of CNOT gates has been presented in the literature. 

The orthogonal gate is an important class of gate, the matrix corresponding to the gate is orthogonal. For example, 
classical reversible logic circuits have a long history |13j and are a necessary subclass whose realization is required for 
any quantum computer to be universal. The matrix elements of them are all real, so they are orthogonal. Utilizing 
C3 the basic property of magic basis, in 2004, Vatan and Williams investigated the synthesis of two-qubit orthogonal 
gate in [7| ■ The result is that the synthesis of the unimodular orthogonal gate requires at most 2 CNOT gates and 
12 elementary one-qubit gates. As for the non-unimodular orthogonal gate, that is its matrix determinant is equal to 
minus one, it requires at most 3 CNOT gates and 12 elementary one-qubit gates [7J]- The number of the one-qubit gates 
required can still be reduced further. Moreover, no articles discuss the synthesis of general orthogonal three-qubit 
quantum gates yet. 

In this work, we devote to investigating the synthesis of two-qubit and three-qubit orthogonal gates. For this 
purpose, we study the Cartan decomposition of matrix for these gates first. Based on the particular decompositions, 
the two kinds of synthesis are obtained. For two-qubit unimodular orthogonal gate, it requires at most 2 CNOT gates 
and 6 one-qubit R y gates, beating an earlier bound of 2 CNOT gates and 12 one-qubit elementary gates. The numbers 
required for one-qubit gate and CNOT gate are all reach the lower bound. For three-qubit unimodular orthogonal 
gate, it can be synthesized by 16 CNOT gates and 36 one-qubit R y and R z gates in the worst case. 
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This paper is organized as follows. The concept of Cartan decomposition and its application in quantum information 
science (QIS) are briefly introduced in Section[nJ Based on a kind of Cartan decomposition of special orthogonal group 
SO (4), we provide an optimal synthesis of general two-qubit orthogonal gate in Section ITTll The decomposition of 
the SO(8) group associated to three-qubit unimodular orthogonal gate is investigated in Section ITVl The synthesis of 
general three-qubit unimodular orthogonal gate is studied in Section [V] It is first time to discuss the synthesis of this 
kind gate. A brief conclusion is made in Section IVII 

II. CARTAN DECOMPOSITION AND ITS APPLICATION IN QIS 

The Cartan decomposition of Lie group Q depends on the decomposition of its Lie algebra. A Cartan decomposition 
of a real semisimple Lie algebra q is the decomposition 

g = iep, (i) 

where p is the orthogonal complement of [ with respect to the Killing form, [ and p satisfy the commutation relations: 

[l,[]C[,[[,p]Cp,[p,p]C[. (2) 

[ is a Lie subalgebra of g. A maximal Abelian subalgebra contained in p is called a Cartan subalgebra of the pair 
(g, I) denoted as a. Then using the relation between Lie group and Lie algebra, every element X of the Lie group G 
can be written as 

X = K X AK 2 , (3) 

where G = e B , K\, K 2 G e 1 and A G e a . 

There are many kinds of Cartan decomposition for semisimple Lie groups. Now the main application in quantum 
information science is the decomposition of SU(2") group for multi-qubit system, i.e. Khaneja-Glaser Decomposition 
(KGD) [J]. Moreover there are some other decompositions, such as Concurrence Canonical Decomposition (CCD) [14, 
[l5j which is a decomposition of SU(2") group too, the Odd-Even Decomposition (OED) [16j, which is a generalization 
of CCD to more general multipartite quantu m sy stem case. Some kinds of Cartan decomposition for a bipartite high 
dimension quantum system were discussed i n I17H19I. These Cartan decompositions have been applied in the synthesis 
and implementation of quantum logic gates [ifl Hll, 120, [2l[ , the entanglement of multipartite quantum systems [14T.I15J. 
etc. But we need to find new suitable algebraic structures of Cartan decomposition to meet the purpose here. 

III. OPTIMAL SYNTHESIS OF GENERAL TWO-QUBIT ORTHOGONAL GATES 

We now consider the decomposition of 4 dimensional special orthogonal group SO (4) associated to the two-qubit 
unimodular orthogonal gate. Difference from that in [7[, the so(4) Lie algebra is constructed as 

so(4) := span{il ® a yi ia y ® I,icr x ® cr y ,i<7 y (g) <j Xl ia z ® <T y , ia y a z }, (4) 

in which each basis vector involves a a y matrix. A kind of Cartan decomposition of so(4) algebra is that 

so(4) = Iep, (5) 

with 

[ := span {il <g> a y , ia v €D /} , (6) 

p := span{ia x <g) a y ,ia y <g) u Xl icr z (g> a y ,ia y <g) a z } . (7) 

where I is a Lie subalgebra and p = [ . Its Cartan subalgebra is 

a := span {ia x ® a y , ia y ® a z } . (8) 

Utilizing the relation between Lie group and Lie algebra, the Cartan decomposition of Lie group SO (4) can be 
obtained. For every element X G SO (4), we have 

X = K X AK 2 , (9) 



where Ki, Ki E SO (2) ® SO (2), and A is a two-qubit operation of the form 

A(a, b) — exp (— i (aa x ® <r y + ba y <g> a z )) , 

where a,b G R. 

The A(a, b) can be represented by the synthesis of elementary gates as 

A = Cl-R v 1 \b)-RV\a)-Cl. 
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FIG. 1: A circuit for computing the general two-qubit unimodular orthogonal gate. 



Here and afterwards Cj denotes the CNOT gate that control on the j'-th qubit and target on the i-th qubit, and 

Ra (a) — exp(—iaa a ) (a € x,y,z) is an elementary one-qubit gate acting on the i-th qubit. Combing Eqs.(9), (10) 
and (11), we can get the synthesis of general two-qubit orthogonal gate as in Fig[Tj it requires at most 2 CNOT gates 
and 6 one-qubit R y gates. As for the non-unimodular orthogonal gate (the determinant is equal to minus one), it 
requires at most 3 CNOT gates and 6 one-qubit R y gates. The 2 CNOT gates is optimal for CNOT gate cost of 
two-qubit orthogonal gate, and it has been proved in [9]. Since a SO (4) matrix has 6 independent parameters, it 
needs at least 6 elementary one-qubit gates to load them. So the synthesis of general two-qubit orthogonal gate here 
is optimal both for CNOT gates and elementary one-qubit gates. 



IV. DECOMPOSITION OF GENERAL THREE-QUBIT UNIMODULAR ORTHOGONAL GATE 

The matrices of any general three-qubit unimodular orthogonal gate are elements of special orthogonal group SO(8). 
We construct the so(8) Lie algebra first. Taking AI type of Cartan decomposition [3| on u(4) and u(2) Lie algebra, 



u(4) = ia {1) © iS {1 \ 



(12) 



u(2)=ia^®zS< 2 \ 
with 

jcr > := span{il <E) <r x , H ® &y, il <£) cr z , i&x 1 1 i&y ® I, i&z ® I}, 

iS^ := span {ia x . y . z ® <J x ,y,z, H} , 

io^ > := span {ia y } , iS^ ' := span{ia x ,ia z ,il}. 
A set of basis for a Lie algebra is given by 28 tensor products of the form 

F := ia (1) <g> S {2) and iS {1) ® a {2) . 
Using the transformation matrix in [22, [2J] : 
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(13) 

(14) 
(15) 
(16) 

(17) 
(18) 



the so(8) Lie algebra can be obtained by M^ ■ span {F} ■ M. So the Lie algebra q := span {F} is isomorphic to so(8), 
and the basis in Eq. (fT5|) can be called as magic basis of so(8) algebra. 
Then we take Cartan decomposition of the Lie algebra g as Eq. ([1]) , with 



[ := span{il ® cr x , y , z ®I,icr x ,y,z ® I ® 7,17 ® Cx.y.z ® &z,i<r x ,y,z ® 7 ® er z }, 

p := span{il ® cr Xi y iZ ® a x ,icr Xt y tZ ® 7 ® a^, £7 ® 7 ® (7 y ,icr x , y , z ® C;z,y,z ® fy}- 
The 1 is isomorphic to so(4) © so(4). The Cartan subalgebra of the pair (g, [) can be chosen as 
o := span{ia x ® a x ® a y , ia y ® <r y ® cr y , ia z ® u z ® a y , il ® I ® cr y }. 



Using the formula 



[A ® £, C ® 7?] = - ({A, C} ® [£, £>] + [A, C] ® {£, £>}) . 



it is easy to verify that the [ and p in Eqs. (fl9l[20|) satisfy the conditions of the Cartan decomposition in Eq.([2]). 
Lie subalgebra [ could be decomposed further 



with 



[=[«©p« 
[W : = span {il ® a x , y , z ® 7, iu x . y . z ® 7 ® 7} , 



pW : = span {?7 ® cr Xi3/i ^ ® cr,,, icr X)1/)jS ® J ® cr z } . 
The [W is isomorphic to so(4). Its Cartan subalgebra can be chosen as 

a^ 1 ' :— span{il ® <7 Z ® a z ,ia z ® I® er 2 } . 



(19) 
(20) 

(21) 
(22) 

(23) 

(24) 
(25) 

(26) 



From the correspondence between Lie group and Lie algebra and the conjugative transformation, we get the Cartan 
decomposition of Lie group SO (8): any element of the group can be decomposed as 



— A/ft st. a^ 1 



(l); 



X = M 1 K 1 AY'K 2 AK 3 A y 2 1 'K 4 M 



(27) 



Here 7Q € SU(2)®SU(2), A and A\ are the Abelian subgroup associated to Cartan subalgebra a and a^ respectively. 



V. SYNTHESIS OF GENERAL THREE-QUBIT UNIMODULAR ORTHOGONAL GATE 

Based on the discussion in Section IIVI the decomposition of general three-qubit orthogonal gate is shown in FigJ^l 
where R = R z (9)R y (tp)R z (^) G SU(2) and 
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FIG. 2: Decomposition of the general three-qubit unimodular orthogonal gate, R — R z R y R z 



A(a, &, c, d) = exp{— i(acr x ® a x ® a v + b<r y ® a y ® a y + ca z ® a z ® <r y + dl ® 7 ® a,,)}, 



(28) 
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FIG. 3: A circuit for computing the magic matrix M. 



A {1 \a, f3) = exp {-i(al <E) cr z ® a z + j3a z ® I ® <r z } . 
The synthesis of transformation matrix .M is given in [7| shown in Fig{3j that is 

*<-<s ■«?»(=) -^G) •«?>(-;). 

The A can be expressed as 

A(a, 6, c, d) = M ■ A(a, b, c) ■ M^ ■ R y 3) (d), 
here 



(29) 
(30) 
(31) 



A(a, b, c) = exp{— i(al <g> a z ® a y — ba z ® a z ® a y + ccr z ® / ® <j y )} 

a-b+c 

b-a+c 

a+b-c 

-a-b-c 



— exp < — i 



(32) 



Since the Cartan subalgebra is commutative, we can break down the synthesis of A(a, b, c) into the following operations: 

Ai(a) = exp{— ial® <j z ® cr y }, (33) 

A 2 (—b)=exp{ibcr z ®a z ®a y }, (34) 

A 3 (c) = exp{— ica z ® I® Cj,}. (35) 
And we have 

Axia) = Cl ■ Rf\a) ■ Cl (36) 

A 2 (-6) = C 3 • Cf • 4 3) (-6) • C 3 • C 3 , (37) 

A 3 (c) = C 1 3 -4 3 )( C )-C 1 3 . (38) 
By putting Eas.(|3"6]). (1371) and ([35]) together, we get 

I(a, 6, c) = Cf ■ 4 3) (a) ■ C 3 ■ R^ (-b) ■ C 3 ■ Rf (c) • C 3 . (39) 
Here the identity C 3 ■ Cf = Cf • Cf is used. Combining Eqs.flHl), (EH and ([3"9"]). we have 

A(a, b, c, d) = C* • if) (J) • C 3 • iif (a) • C 3 • i?< 3 > (-6) ■ 

Cl-R y ^{c)-Cl-R^(-\)-Cl-Rf(d). (40) 

and its circuit shown in Fig|4j Since R z gates commute with the control qubit of the CNOT gate, here the R z gates 
in M. and M^ are canceled. The synthesis of A^ is shown in FigJSl that is 



A«(a,/3) = C 3 • 4 3 )(/3) • Cf • C 3 • i?i 3 >(a) • C 3 



(41) 



Putting all these pieces together, we get that 16 CNOT gates and 36 one-qubit R y and R z gates at most are 
sufficient to synthesize general three-qubit orthogonal gate. For the same reason mentioned above, the two R z gates 
acting on same qubit neighbored A^> have been combined to one. Here, each M. requires 1 CNOT gates and 3 
one-qubit R y and R z gates, A requires 6 CNOT gates and 6 one-qubit R y gates, each A^ requires 4 CNOT gates 
and 2 one-qubit R y and R z gates, and 8 R gates require 20 one-qubit R y and R z gates. 
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FIG. 4: A circuit for computing the unitary operation A. 




FIG. 5: A circuit for computing the unitary operation A- 1 ' . 



VI. CONCLUSIONS 



Based on the decomposition of matrices, the synthesis of two-qubit and three-qubit orthogonal gates is investigated. 
For two-qubit orthogonal gate, we get optimal result, which requires at most 2 CNOT gates and 6 one-qubit R y gates, 
beating an earlier bound of 2 CNOT gates and 12 one-qubit gates. For the three-qubit unimodular orthogonal gate, 
it requires 16 CNOT gates and 36 one-qubit gates from the family {R y ,R z } in the worst case. There are abundant 
algebraic structures for matrix decomposition of three-qubit orthogonal gate. We have many ways to investigate the 
synthesis of general three-qubit orthogonal gate. The result given here is the best one we have got, although we can 
not affirm that is optimal yet. The synthesis of general three-qubit gate has been studied in some literatures |H>ll2l |. 
the orthogonal gate is an important class of gate of them. So the work here is essentially on the "small circuit" 
issue of three-qubit gates, which is first investigated in this paper. Different from two-qubit gate, how to get optimal 
quantum circuit for general three-qubit gate has not been well solved and is worthy studying further. 
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